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The Trojan-Horse method is an indirect approach to determine the energy dependence of S- 
factors of astrophysically relevant two-body reactions. This is accomphshed by studying closely 
related three-body reactions under quasi-free scattering conditions. The basic theory of the Trojan- 
Horse method is developed starting from a post-form distorted wave Born approximation of the 
fSJ , T-matrix element. In the surface approximation the cross section of the three-body reaction can 

. be related to the S-matrix elements of the two-body reaction. The essential feature of the Trojan- 

' Horse method is the effective suppression of the Coulomb barrier at low energies for the astrophysical 

^N) ^ reaction leading to finite cross sections at the threshold of the two-body reaction. In a modified 

plane wave approximation the relation between the two-body and three-body cross sections becomes 
very transparent. The appearing Trojan-Horse integrals are studied in detail. 
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Q> I. INTRODUCTION 

\D ■ 

] Nuclear reaction rates are an indispensable ingredient of many astrophysical models. They have to be known with 
' sufficient accuracy in order to understand quantitatively the evolution of the universe, stars and other objects in the 
cosmos 1^, ^, Q. Ideally, reaction cross sections are directly measured in the laboratory. But, with few exceptions and 
despite many experimental efforts, the relevant low energy range cannot be reached in direct experiments ^, 0- 
Cross sections for reactions with charged particles rapidly become very small with decreasing energy of the colliding 
, nuclei due to the repulsive Coulomb interaction. Extrapolations of the cross section a{E) to low energies from results 
' I ' at higher energies accessible to experiments are often needed. This is accomplished with the help of the astrophysical 
O ' S-factor 

Ch ■ S{E) = a{E) E exp(27r?7) (1) 



where E is the cm. energy and r] = Z\Z2e^ j {hv) is the Sommerfeld parameter which depends on the charge numbers 
Zi , Z2 of the colliding nuclei and their relative velocity v in the entrance channel. The S-factor shows a much weaker 
energy dependence than the cross section a-{E) because the main effect of the penetrability through the Coulomb 
barrier is compensated by the increase of the exponential factor. The extrapolation process introduces uncertainties 
and important contributions to the cross sections, like resonances, can be missed. Additionally, direct laboratory 
experiments are affected by electron screening, which effectively reduces the Coulomb barrier between the nuclei and 
enhances the measured laboratory cross section [|[ ^, |l^ . A correction has to be applied to obtain the cross section 
for bare nuclei. This effect does not seem to be completely understood yet and independent information on low 
energy cross sections is valuable in order to develop a quantitative description of electron screening. In astrophysical 
applications one has, in addition, to account for the screening under stellar conditions. 

In recent years several indirect methods have been developed to extract cross sections relevant to astrophysics from 
other types of experiments. In these alternative approaches the astrophysical relevant two-body reaction is generally 
replaced by a suitably chosen three-body reaction. The relation between the reactions is established with the help 
of nuclear reaction theories. Without doubt, this process will introduce some uncertainties, but valuable information 
on the astrophysical reaction can be obtained. Also, the errors are independent from that of the direct approach. Of 
course, the indirect methods have to be validated by studying well known reactions before firm conclusions can be 
drawn from indirect experiments in cases where direct measurements are not feasible; see also the minireview pT| |. 

The Coulomb dissociation method has become a valuable tool for extracting low-energy cross sections of radiative 
capture reactions 0(6,7)0 by studying the inverse process of photo dissociation 0(7,6)0 Instead of using real 

photons, the Coulomb field of a highly charged target X acts as a source of virtual photons during a scattering 
process which leads to the breakup reaction X{c, ab)X with three particles in the final state. Due to the high flux of 
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virtual photons, cross sections at the small relevant energies in the two-body system are strongly enhanced as compared 
to the capture reaction. Another approach is the ANC method that tries to extract the asymptotic normalization 
coefficient of a nuclear ground state wave function by studying various combinations of transfer reactions involving this 
nucleus at low energies jlj, 0. The coefficient can be used to determine theoretically the astrophysical S-factor for 
radiative capture reactions at zero energy. However, these indirect approaches are limited to astrophysically relevant 
reactions where a photon is involved. 

The observation of a similarity between cross sections for two-body and closely related three-body reactions under 
certain kinematical conditions led to the introduction of the Trojan-Horse method (THM) see also [0 p^ . 
The aim of the THM is to extract the cross section of an astrophysically relevant two-body reaction 

A + x-^C + c (2) 

from a suitably chosen reaction 

A + a-. C + c + h (3) 

with three particles in the final state assuming that the Trojan Horse a is composed predominantly of clusters x 
and b. The kinematical conditions are chosen such that the momentum transfer to the nucleus b is small during the 
reaction. Therefore b can be considered as a spectator to the reaction of A and x. This process is often referred to 
as a quasi-free scattering. In a selected part of the available three-body phase space it is known to dominate over 
other reaction mechanism like sequential breakup processes. In the past quasi-free scattering has been used to extract 
information on momentum distributions of the nucleus a, i.e. the ground state wave function in momentum space, 
employing a plane- wave impulse approximation (PWIA) in the theoretical description [ p^ . In this approach the cross 
section of reaction (|^) factorizes into a kinematical factor, the ground state momentum distribution of nucleus a and 
an off-shell cross section of reaction (^) that is assumed to be known. On the other hand the two-body cross section 
can be extracted from the cross section of reaction (|^) if the momentum distribution of the Trojan Horse a is known 
with sufficient accuracy and a relation between the off-shell and on-shell two-body cross sections can be established. 
The selection of different spectators b and thus Trojan Horses a allows additional checks of the underlying assumptions 
of the method. 

Unlike the Coulomb dissociation method and the ANC method which are limited to radiative processes, the THM 
can be applied to reactions where no photon is involved. The essential feature of the THM is the actual suppression 
of the Coulomb barrier in the cross section of the two-body reaction. The cross section of the three-body reaction 
is not reduced when the cm. energy in the A -\- x system approaches zero as in reaction (H). The energy in the 
entrance channel of (|^) can be around or above the Coulomb barrier and effects from electron screening are negligible. 
Nevertheless, very small energies in the reaction (||) can be reached. 

The feasibility of the THM was studied in several experiments involving various reactions during the last years 
pO| , [2l| , p2| , [2^ , p^ , [25I , p6| , [2^ , pst . In earlier evaluations of the experiments some simplifying assumptions were 
made in the theoretical description. The off-shell two-body cross section as appearing in the PWIA was considered 
as the bare nuclear cross section. It was converted to the on-shell two-body cross section by correcting for the 
Coulomb penetration in a heuristic approximation. Basic considerations for the theoretical description of the THM 
were developed in ||2^ but the deduced relation between the cross sections for reactions (|^) and (H) were not directly 
applicable to the experiments. 

In this paper the theory of the THM is developed in certain approximations that allow to establish a simple 
connection between the cross sections of the three-body and two-body reactions. In Section |l| the reaction theory 
is formulated and the relation of the T-matrix element of the three-body reaction with the S-matrix elements of the 
two-body reaction is found. In connection with a plane wave approximation, fundamental TH integrals appear that 



are discussed in Section III. In the following two sections expressions for scattering amplitudes and cross sections 
are derived for spinless particles and then for the general case of particles with spin. The TH Coulomb scattering 
amplitude that is relevant in the indirect study of elastic scattering processes is treated in Section Kinematical 



conditions, the energy dependence of cross sections and applications ot the THM are discussed in Section VII. A 
summary and an outlook are presented in the last Section. Details of the analytical calculations are given in the 
Appendices. The evaluation of actual Trojan-Horse experiments with the present theory is beyond the scope of this 
paper und is subject to a detailed treatment in separate studies. 



II. REACTION THEORY 



In the three-body reaction (||) the particle x is transferred from the Trojan Horse a to the nucleus A leading to the 
C -f c final state with the spectator b of the Trojan Horse remaining, see Fig. |l]. This reaction can be considered as 



FIG. 1: Transfer of particle x in the Trojan-Horse reaction (g). 



a transfer to the continuum, where the inelastic reaction A + x ^ C + c can happen in addition to the elastic A + x 
scattering. It is customary to describe such a direct reaction, i.e. A + a— > C + c-l-6, with the help of a distorted wave 
Born approximation. It is not required in the theoretical approach that A^- x ^ C -\- c] elastic scattering processes in 
the two-body system can be treated in the formalism, too. Effects from the antisymmetrization of the wave functions 
will be neglected in the present treatment; they are expected to be small. 



A. Coordinate Systems and Cross Sections 

In a three-body system various sets of Jacobi coordinates are used to specify the positions of the particles. In the 
theoretical description we encounter the sets 

^ ^ ^ ^ ^ ^ ^ m^r^+nibn , 

r-xb ^r^-n , TAa = - Ta = TA ; (4) 

rux + mb 

in the initial partition and 

^ ^ ^ ^ ^ ^ mcfc+mcrc ^ 

rcc = rc -Tc, TBb = rs - n = ; n (5) 

mc + rtic 

in the final partition. The symbol B denotes the C + c — A + x system. The coordinate vectors are shown in Fig. ||. 
The corresponding relative momenta and wave vectors are given by 

= nk., = HhE,^ (6) 

for nuclei i and j with masses to^, mj and momenta Pi,Pj in the laboratory system, respectively. With the help of the 
kinetic energies 

E 

where the reduced masses 

- (8) 
trii + rrij 

appear, energy conservation in the two-body reaction can be expressed as 

Eax = Ecc - Q2 (9) 

with the Q-value 

Q2 = {ruA + mx - mc - mc)c^ (10) 
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and similarly 



EAa — Ecc + EBb — Qa 



Ea 



E 



Bb 



Q2-Q3 



with 



(11) 



(12) 



in case of the three-body reaction. 

The general cross section for reaction (||) with three particles in the final state depends on the choice of the 
independent variables and the reference system. In the cm. system the energy Eqc and the directions of the Jacobi 
momenta pBb and pcc completely specify the kinematical conditions for a given projectile energy. Then the differential 
cross section takes the form 



dEccdVLccdVLst 



Tfi{kcc, ksb] kAa] 



with the cm. kinematical factor 



^J-Aa^J-BbfJ-Cc PBbPCc 



PAa 



(13) 



(14) 



and the T-matrix element Tfi. In an actual experiment, the particles C and c are usually detected in the laboratory 
system and it is more convenient to use 



dEcdVLcdn, 



Tfi{kcc, ksb] kAa) 



with the kinematical factor 



Klah 



P-Aafnc PCPc 





_ PCc\ 


Pc 


\pBb 


mc J 


Pc_ 



(15) 



(16) 



The laboratory differential cross section depends on the scattering angles of the nuclei C and c and the energy Ec- 
Again, these quantities, together with the beam energy, specify the kinematical conditions completely. The nuclei A 
and a can be projectile and target or vice versa. In case of particles with spin appropriate averages over initial states 
and sums over final states have to be considered. 



B. Approximations of the T-Matrix Element 

The T-matrix element Tfi in ( p^ ) and (|l5|) contains all the essential information relevant to the scattering process. 
It has to be calculated in a suitable approximation that allows to find the connection to the cross section of the 
astrophysical relevant reaction (p|) . 
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The Trojan-Horse reaction ^ has the form of a usual two-body reaction 

A + a-^B + b (17) 

if the C + c system is considered as an excited state of the compound system B. The exact T- matrix element for this 
reaction is given in the post-form description by 

Tfi = (exp (^iksb ■ rsbj (l)B(t>b\VBb\^A^ (kAa^rAa)) (18) 

with the exact scatttering wave function "^^y^J {kAa,'rAa) in the initial state and the potential Vsh between B and 
b in the final state ||30|] . The wavefunction 0^ of a nucleus i depends only on internal coordinates which are not 

given explicitly. The relative motion of B and b is described by a plane wave with momentum hhsb- Applying the 
Gell-Mann Goldberger transformation the T-matrix element (^8|) assumes the form 

Tf^ = {XBbikBb,rBb)Mb\VBt - UBb\'i''-Aj{kAa,rAa)) (19) 

which is still an exact relation. For a detailed derivation see Here, X^m (^BbT^Bb) is a distorted wave for the B — b 
relative motion generated by a suitably chosen optical potential UBb that only depends on rBb but not on internal 
coordinates in contrast to Veb- In the THM the wave function (j)B does not describe a bound state but a complete 
scattering state 

c^B^-^^chkcc.rcc) (20) 

since the system _B = C + c is in the continuum. 

In the post-form distorted-wave Born approximation (DWBA) the T-matrix element (^9|) is replaced by 

Tf^^'^ikccMb-Ma) - {XBhkBb.rBb)-^^chhcc.rCc)MVxb\XAl{kAa,rAa)<l>A<l>a^ (21) 

where the exact wave function 'I'^''(fcAa, ?Via) is replaced by a distorted wave X^AaO^AairAa) for the A — a relative 
motion. Additionally, the potential Veb ^ UBb is approximated by V^b |2^, This is the usual starting point of 
T-matrix calculations for direct transfer reactions. 

In the next step, the so-called surface approximation is applied which is essential to the THM. For small distances 
between the colliding nuclei the optical potentials are usually strongly absorptive and only reactions at the surface of 
the nuclei contribute significantly to the matrix element. Therefore the full scattering wave function "^[yj (kccT^cc) 
can be replaced by its asymptotic form for radii larger than a suitably chosen cutoff radius R that is larger than the 
range of the nuclear potential. It is typically in the range of the sum of the radii of the two colliding nuclei. The 
interior part of ^P^J (fccc ^^s'c) is set to zero. The validity of the surface approximation was checked in it was 
found to be quite good for the (d,p) reaction at Ed = 26 MeV. The asymptotic form of the scattering wave function 
for Va > R, where a = Ax, Cc, . . . is a possible partition of the system B, is given by 

*gasym(fccc, rcc) = ^ E E (^^ )r,:;, (fccc)0a (22) 



kcc , 

a In 



with radial wave fuctions 



— ^; 

where the Coulomb wave functions 13 



£!i^\a) = /— S'-^(j^uf \rja]kara)-5aCcUl \r]a\kara) , (23) 
t\c.) = il+^*{a) (24) 



[±)(,y;fcr) = e^""'-'^^[Gi{T];kr)±iFi{T];kr)] (25) 

/ It: 
—> exp zLi Ikr ~ T] ln(2A:r) — — 



appear. The Sommerfeld parameter 



V^J = (26) 
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depends on the charge numbers Zi, Zj of nuclei i, j and their relative velocity vi^ = Pij/^ij in the partition a = ij. The 
Coulomb phase shifts are given by (Ji{ri) = argr(l + I + irj). The S-matrix elements Sl^(~.^ in the radial wave functions 
completely describe the two-body scattering process. They are on-shell quantities and the momenta ka are derived 
from energy-conservation in the two-body reaction. E.g., in the a ~ Ax partition we have kAx — V^Jj^AxEax / with 
the energy Eax from relation (H). 

The potential Vxb appearing in the T-matrix element ( pi]) describes the interaction between x and b in the Trojan 
Horse a. Assuming a simple cluster picture of a and neglecting contributions from excited states of x and b or other 
partitions in the ground state, the momentum amplitude W{q) of the product 

f d^Q 

Vxbirxb)4'a{rxb) ^ J j^;^W{q)exp{iq-rxb)'f'x(t)b (27) 
can be introduced and the T-matrix element assumes the form 

Tj/'ikcc, kBb; kAa) = j ^ W{q) (28) 
iXmi^Bb, ^S6)*c"c!asym(^Cc, ^Y^c)! exp (iq- f^b) X^Aai^Aa, rAa)<pA<Px) ■ 

The assumption W = const, corresponds to the zero-range approximation which is frequently employed in the calcu- 
lation of DWBA T-matrix elements. In this case, the integration over q leads to a (5-function in the variable f^b and 
the actual dimension of the appearing integral is reduced from six to three. However, in general the full momentum 
dependence of the amplitude W has to be considered. With the help of the Schrodinger equation the momentum 
amphtude W is related by 

Wiq)^-(^ea + ^^<faiq) (29) 

to the ground state momentum wave function 

*a(g) = (exp {iq- rxb) 4'x(t)b\4'a{rxb)) (30) 

of the nucleus a. The energy £a = {mx + rub — ma)c^ = Q2 — Q3 > is the binding energy of a with respect to the 
X + b threshold. 

The integration over the internal coordinates in the matrix element ( p8|) selects the a = Ax partition of the 



asymptotic wave function (22) and the T-matrix element in the TH approximation can be written as 



Tj/'ikcckBb^kAa) - :p^\/^E(2^ + l) (31) 
■' 2lkcc V VAx ^ 

SaxCc Ui^^ [ksb, kcc, kAa) - SaxCc C// ^ [ksb, kcc, kAa) 
It resembles a scattering amplitude of the reaction C + c ^ A + x except for the functions 

U^^\kBb,kcc,kAa) ^P^^T.YUkcc) f 70 Wiq) (32) 

iX^Bbi^Bb, TBb)^^^^ —Ui^\'nAx;kAxrAx)Ylni{fAx)\GX-P {iq ■ rxb) x'ita i^Aa, TAa)) 

rAx 

that describe the angular distribution and the momentum dependence due to the presence of the spectator h in the 
reaction (^. The unit step function 9 accounts for the surface approximation by eliminating the contributions at 
small radii. With equation ( |3l| ) the relation between the cross section (|l5|) of the three-body reaction (||) and the 
S-matrix elements of the two-body reaction is directly established. It is possible to extract the energy dependence of 
the S-matrix elements from experimental three-body cross sections, at least in principle. For an actual application 
of the TH method, it is convenient to introduce additional approximations that lead to a formulation with a direct 
relation between the cross sections of the two-body and three-body reactions similar to the PWIA. 

Even with the surface approximation, the remaining matrix element involves a six-dimensional integration in the 
Jacobi-coordinates which is an extensive computational task. Various approximation can be introduced to simplify 
the calculation. It is useful to expand the distorted wave in the initial state in a Taylor series 

X^a' i^Aa, TAa) = + i^^Aa " TAx) • V + . . . ) X^aI {^Aa^rAx) (33) 

« exp i {rAa - rAx) ■ K X^Aa i^Aa,rAx) 
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where the wave vector K replaces the derivative — iV with respect to the spatial coordinates. In the so called local 
momentum approximation ||3^, the modulus of K is determined by the actual kinetic energy of the A — a relative 
motion at a certain distance Rau- The direction of K is assumed to be same as the asymptotic momentum kAa- If 
XaJ is ^ plane wave and K = kAa the relation (|3^) is obviously exact. Considering the relations 

rxb = rsb - arAx , rAa = fifsb + (1 - OiP)rAx (34) 
of the Jacobi coordinates with factors 

a = ■ , /3 = ■ , (35) 

rriA + mx mf, + 

the matrix element in equation (|2^) factorizes and the T-matrix element in the TH approximation assumes the form 

fq_ 

(2nf 



Tj^"{kcc, kBb\ kAa) = / 7^ W(g, kBb) M{q, kcc, kAa) (36) 



with the generalized momentum amplitude 

Wiq, kBb) = Wicf) {xm (kBb, rBb)\ cxp (^Q • VBb) ) (37) 

and the matrix element 

M{q, kcc, kAa) = (*Cc!asym(^Cc, fcc)\exp (^-iaQ ■ TAx^ X^Aai^Aa-, rAx)(l> A(l>x) (38) 

with Q ^ q + I3K. The integration over the internal coordinates leads to the expression 

M{q,kcc,kAa)^T^.[^y^{2l + l) (39) 

llKCc V '^Ax J 



SaxCc^I iq,kcc,kAa) - SaxCcI^i {q,kcc,kAa) 



(±),.r^ ? N 4^(-»)' [ 0{rAx-B)(±). 



with the integrals 

l^l^'{qikcc,kAa) = TKTT-T^ I '^'^'^Ax"-^-^ —ufi-qAx^kAxrAx) (40) 

[IL + 1) J rAx 

exp {-iaQ ■ rAx'j XAi^i^Aa, tax) ^ Yi*^{fAx)Yim{kcc) 

m 

The functions C//"^^ in the TH T-matrix element (|l|) assume the form 

(27r)3 



Ul^\kBb, kcc, kAa) - / W(g, kBb) U\^\q, kcc, kAa) ■ (41) 



The dependence of L^/^' on kcc or equivalcntly kAx together with the energy dependence of the two-body S-matrix 
elements SK^q^ leads to a finite cross section of the three-body reaction ^ even when the threshold of the two-body 
reaction (||) is reached. This is the essential feature of the TH method. However this is not readily seen from the 
general expression (p2). A simplified formulation allows to study the energy dependence more explicitly. 



C. Plane Wave Approximations 

If Xsb is assumed to be a plane wave, the matrix element in the generalized momentum distribution ( |37| ) leads to 
a (5-function in the variable q 

>V(g, kBb) = W{q) (27r)3 ,5 {q- Qeb) (42) 
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with 

Qsb = ksb - PK . (43) 



Then, the q-intcgration is trivial; the T-matrix clement factorizes explicitly into a momentum amplitude and the 

matrix element ( ^9| ) where the integrals ( ^ ) are evaluated for q = Qst- This plane- wave approximation for xhb 
justified when the energy of the Bb relative motion is large. When the spectator 6 is a neutron, it is also often possible 
to replace Xbi! ^ plane wave. On the other hand, when x^Bb replaced by a pure Coulomb scattering wave, e.g. 
if the Bb relative energy is small, the matrix element in (|3^) can be calculated explicitly (see Appendix |b|). However, 
in applications of the THM the cm. energy in the Bb system is generally large because of the large projectile energy 
and Q- value of the reaction and the latter case is not relevant. 

Independently from the treatment of the wave function x^sb ^ ^ plane wave approximation can be introduced for 
XAa ^^'^ matrix element (^). In this case one finds 

M{q, kcc, kAa) = (*Cc!asym(^Cc, rCc)| exp (iQAa ' TAx^ 4>a4>x) (44) 



with 



The quantities (Efl) reduce to 



QAa = kAa ^ a(3K - aq . (45) 



Ul^\q,kcc,kAa) = T~k ^liQAa ■ kcc) Ji^\R,'nAx,kAx,QAa) (46) 

with Legendre polynomials Pi and the dimensionless Trojan-Horse integrals 

oo 

4^\R,VAx,kAx,QAa) = kAxQAaJdrAxUr\vAx;kAxrAx)rAxMQAarAx). (47) 

R 

The Coulomb wave functions as defined in (|2^) and the regular spherical Bessel function ji of order I appear in 
the radial integral. It is convenient to introduce the decomposition 

4^\R. VAX, kAx, QAa) - e^^^'f"^-) [Ji^HR, VAx, kAx, QAa) ± ^/^^(i?, TjAx, kAx, QAa)] (48) 

of the TH integrals into contributions with the regular and irregular Coulomb wave functions. The integrals are 



discussed in detail in Section HI 



In the following only the plane-wave approximation for both Xau ^^'^ Xhb ^^^^ considered. In this case the 
T-matrix element has the simple product form 

Tj,"{kcc,kBb;kAa) = W{QBb){^'hl.,yAc,rcc)\(^^v[^QAa-rAx)Mx) . (49) 

This approximation already contains the essential ingredients to see the principles of the TH method and the connec- 
tion to the PWIA becomes clear. Generalizations to a more general treatment with distorted waves are obvious from 
the above. The argument of the amplitude W is the momentum 

Qsb = kBb ■ kAa (50) 

nib + nix 

assuming that K — kAa- Neglecting the Fermi motion of 6 inside the Trojan Horse, mbkAal {'TUb-'r'mx) is the momentum 
of b relative to A in the initial state and kBb is the momentum of b relative to B = C -I- c in the final state. Thus 
—QBb corresponds to the momentum transfer to the spectator b. The amplitude W describes the distribution of the 
transfered momentum due to the Fermi motion. Similarly, the momentum 

QAa = kAa ^Bb (51) 

rriA + nix 
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in the argument of the plane wave can be considered as the (negative) momentum transfer to nucleus A (independent 
of the choice of K) by the particle x. In the case of a infinitely heavy nucleus A, equations (^0|) and ( |5T| ) reduce to 

Qsh = ] ka kfj — kx ; ka , QAa — kiy ka — kx (52) 

mfc + nix ITT-b + ITT-x 

(cf. Fig. |l|) and the interpretation becomes simpler. The main task is to calculate the TH integrals (^) and to find 
the explicit relation of the three-body cross section to the two-body cross section. 

The plane wave approximations for Xai^ E^^d XBb seem to be crude at first sight but the TH matrix element ( |49| ) 

still contains the asymptotically correct wave function ^'irym(C'c, fccc) for the C + c system and thus the complete 
information on the two-body scattering process. This is in clear contrast to the PWIA where the effect of the 
Coulomb barrier in A + x system on the energy dependence is not obvious. The absorptive feature of the optical 
potentials is taken into account by the surface approximation. Additionally, the distorted wave XaJ describes a 
scattering state with a much higher (and constant) momentum as compared to the A + x system and only a small 
part of the three-body phase space is of interest in the reaction. The absolute cross section for the three-body process 
calculated in the plane wave approximation might be different from the actual value but the energy dependence 
with respect to the two-body reaction is expected to be treated correctly. 



III. TROJAN-HORSE INTEGRALS 



The main difference between the usual expressions for the scattering amplitude and the corresponding T-matrix 
element (^l|) in the TH method is the appearence of the functions U^^^ . They decribe the off-shell behaviour of the 
two-body scattering amplitude and lead to the effective removal of the Coulomb barrier in the quasi-free scattering 
process. This effect can only be understood if the energy dependence of these functions is known. In the plane-wave 
approximation the discussion becomes very transparent. The functions U^^'^ factorize into a momentum amplitude, 
Legendre polynomials and the TH integrals ( ^7| ) which arc not difficult to study. 



A. General form of Trojan-Horse integrals 



The basic Trojan-Horse integral in the plane wave approximation is given by 



4"\R,7^,k,Q) 



dr Hi{ri; kr) zi{Qr) 



(53) 



where Hi is a Coulomb wave function Fi or G; and zi{x) — xji{x) is a Riccati-Bessel function. Then the integrals j/^'' 
are obtained from equation (|8|). The TH integrals do not converge in the usual sense, since Hi and zi oscillate with 
constant amplitude if r goes to infinity. Convergence is achieved only in the distributional sense after an integration 
over Q with a suitable test function. The problem is caused by the fact that the matrix element in eq. ( |2^ ) contains 
the overlap of continuum wave functions. A similar problem occurs in the case of Bremsstrahlung matrix elements 
that have to be regularized. 

However, the integral ( p3|) can be transformed into a form that allows a numerical calculation. With the differential 
equation for the Coulomb wave functions 







k^ dr^ 



Hi{kr) 



the TH integral (M) becomes 



2?7 
kr 



Hi{kr) 



(54) 



OC 

4"\R,'n,k,Q) = k j dr 



2ri 1(1 + I) _ 1 d^ 
kr k^r^ k^ dr'^ 



Hi{ri\ kr) zi{Qr) 



(55) 



After partial integration and with help of the differential equation for zi 

l{l + l) 



1 d2 
0=Q2^-'(Q-) 



1 



ziiQr) 



(56) 
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the expression 



B ~ Q2 



(57) 



2r)kll"\R, ry, k, Q) + kH[{f^- kR)zi{QR) - QHi{f]; kR)z'i{QR) 



with the converging integral 



/p)(i?, 77, fc, Q) - j dr HiiTj; kr) r"! zi{Qr) 



(58) 



is obtained. The prime denotes the differentiation with respect to the argument. We disregard the contribution from 
the upper bound (infinity). This corresponds to the regularization procedure mentioned above. Relation (^7|) c annot 
be used if fc = Q but this case will never occur since k < Q from kinematical considerations (see Section VII A ). The 
integral (|5^) can be written as 



IX 

4"\R,r],k,Q) = kQMi-\H,r],k,Q)~ J drHiirj;kr)r-^ ziiQr 



(59) 



with the integral 



M,;,^JiJ,,7,fc,Q) = ^J drHi,{j^-kT)r-^zi,{Qr) 



(60) 



because is still finite in the limit R ^ 0. The functions M^^^^ are special cases of the radial matrix elements in 
the quantum mechanical theory of Coulomb excitation ]3^ . A more general form also appears in the general theory 
of transfer reaction to the continuum |Q. In the THM only monopole matrix elements (A = 1) are needed. They 
can be calculated explicitly for both the regular and irregular Coulomb wave functions (see Appendix^ for details). 
The TH integrals with the regular Coulomb functions are given by 



Moo\F,Tj,k,Q) 
MriHF,rj,k,Q) 



Coiv) 



sinC , 



2r]kQ 
irikQ{l + r]^) 



k Q 



sin C ~ 277 cos C 



(61) 
(62) 



for ^ = 0. 1. The argument of the trigonometric functions is = rj In zi with zi — (Q + k)/{Q — k). The constants 
Ci (77) are recursively defined by [p3 



2Tn] 



(p , 2\ 2 



,exp(27rr7) - ly ' " 1(21 + 1) 

The case with the irregular Coulomb wave functions is slightly more intricate and leads to 



M^o^{G,T],k,Q) 
M^,'{G,ri,k,Q) 



2'qkQ ^ 2fcQCo(?7) 
3Ci(77) 



Ar]kQ{l + rf) 

TT 

^12fcQCi(77) 



^ + g ) (l-cosC)-2r/sinC 



(63) 

(64) 
(65) 



^ + -]{l-S,) + 2r^S, 



with the series 
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Explicit expressions for larger I are obtained from the recursion relation {I > 0) 

^ + f ) ' = 1^ + + + l + (67) 

which is valid for Hi — Fi and Hi — Gi. It can be derived in a similar manner as in Coulomb excitation theory [ p7[ . 
However, for a numerical calculation the question of stability of (|67| ) arises. In the case of Hi — Fi only the backward 
recurrence is stable. Starting with, e.g., M£l{H, rj, k,Q) = and M£^-^j^_^{H, ry, k, Q) — 1 for large I a, backward 
recursion gives integrals for I = 0, 1, 2, . . . after normalizing to the known value (p^). In the case of Hi = Gi the 
upward recurrence can be used with the known starting values Mqq^{G, rj, k, Q) and M^^{G, rj, k, Q). The remaining 
integral from to i? in equation (M) is easily performed numerically. 



When T] becomes large (and k small) it is not favorable to calculate the integral {R, ry, fc, Q) from equation (|5^ 
The irregular Coulomb wave functions becomes very large for small radii r and a difference of large numbers has to be 
evaluated with loss of accuracy. In this case it is more convenient to use equation (|5^) directly because Gi becomes 
small with increasing r > R very rapidly. 

B. Limiting cases 



In the case of i? ^ the surface contributions in eq. (|57| ) vanish only if Hi = Fi. This can be seen from the 
approximations 



Piiv,^) - ^^^^^^^V^l2i+i{2V2^) « Ci{v)x'+' (68) 



and zi{x) — ^ a;'+^/(2Z + 1)!! for x = kR 0. In this case the expression for TH integral with the regular Coulomb 
wave function reduces to the simple form 

jP{0, 77, fc, Q) ^ ^!^M[^\Fi,rj, k, Q) . (69) 

Using the approximation Q 



for small x one obtains for Hi = Gi the result 



2TjkHln\Gur],k,Q) 



(2; + l)!!G(ry) 



(71) 



with a remaining contribution from the surface terms. 

In the special case where the transfered particle x is a neutron, the Coulomb functions ul^\x) reduce to Hankel 
functions x[—yi{x) ± iji{x)\ = -iiixh^j^\x) with the spherical Neumann and Bessel functions yi and j;, respectively 
P^ . There is no contribution from the term with the integral I^^^ and the TH integral simplifies to 

j\^\R,Q,k,Q) - ^^-^ [kh\^l{kR)ji{QR) - Qh^^\kR)]i^,{QR)] (72) 

so that Ji~^\R,0,k,Q) is ifc^Q times the J; integral in the stripping enhancement factors Fi of Ref. For the 

neutral particle, there is, for / 7^ 0, still a barrier, the angular momentum barrier with its l{l + l) jr^ radial dependence. 

The enhancement is due to the Hankel functions [x) with a x^'-'^^-' behaviour for a; ^ 0. The mathematical 
treatment of the Coulomb functions is much more involved. 

C. Energy dependence of Trojan-Horse integrals 



By a simple change of variables it is seen that the TH integrals (|53| ) depend only on three independent parameters. 
With 

C = ^ and xo = QR (73) 
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we can write 



4"\R,V,k,Q) ^ 4^J^{xo,V,0 = ^ / dxHi{r,-ix) zi{x) 



(74) 



with reduced TH integrals that are functions of xq, 77, and ^. In order to investigate the behaviour of J( for A; — > 
or ^ ^ it is useful to study the case xq = 0, i.e. i? = 0, first. From the explicit expressions for the integrals M^^{Fi) 
the approximations 



jf)(0,77,fc,Q) = J^"^^,(0,77,O ^ -Co(?7)^"sin(2770, 



4^)(0,ry,fc,Q) = j|;^^, (0,77,0 



sin(2?70 



cos(27/0 



(75) 
(76) 



for the integrals with the regular Coulomb functions and Fi are found. Because the recursion relation (|6^) for 
MjJ^ reduces to the recursion relation of the Riccati-Bessel functions we have the approximation 



(77) 



for all I in the limit fc ^ 0. The argument 277^ = 2rjk/Q is independent of k for constant Q and the entire k dependence 
is determined by Co{ri)S^'^ independent of I. 

Figure Ijshows the dependence of the scaled reduced TH integrals jl^^^iO, rj^ ^^vivO ^i^^ / = 0, 1 on ^ = k/Q 

for various values of r/^ = Z aZ ijlax / {fi?' Q) which is a constant for given Q. The scaled TH integrals become constant 
for small ^ as expected. The rapid decrease of the unsealed TH integrals with decreasing ^ for constant 77^ is obvious. 
The angular momentum barrier leads to a smaller TH integral with larger I for the same parameters ^ and 77O This 
effect is more pronounced for small 7/^ than for large 7/^ where the Coulomb barrier dominates the ^ dependence. The 
approximations ( |75| ) and (^6|) agree very well with the exact results for the TH integrals with Z = 0, 1 for large 77^ and 
not too large The agreement is less satisfactory for small 77^ and larger I as one may expect, since then the angular 
momentum barrier is getting more important as compared to the Coulomb barrier. 

The TH integrals with the irregular Coulomb wave functions reduce in the limit fc to the form 



with functions 



jf)(0,7;,A;,Q) = j(;;)^(0,77,e) 



jf)(0,7;,A;,Q) = jg^^(0,77,O 



Co(r/) 



[1 - 277$ /(2?7e)] 



Co(^) L'yC 



-/(277e)- 277$ 5(2770 



(78) 
(79) 



6 {IT 

f{x)= / dt = Ci(a;)sina;- [^Si(a;) - - 



cos a; , 



g{x) ^ dt 



te- 



t2 + l 



-Ci(a;) cosa; — Si(a;) — — 



smx 



(80) 
(81) 



that can be expressed in terms of the sine and cosine integrals Si(2:) and Ci(a;) [Q. The expressions inside the 
brackets in eqs. ( [zs] ) and (|79| ) are constants and the fc-dependence of the TH integrals with the irregular Coulomb 
wave functions is determined by £,/Co(j])- 

In Figure ^ the dependence of scaled reduced TH integrals ^^/redC^i exp(— 7777) with / = 0, 1 on ^ is shown 

for the same range of the parameter 77^ as in Fig. ||. For small ^ the scaled TH integrals again become constant. 
The unsealed TH integrals increase dramatically for decreasing ^ for constant 77^ and for increasing 77^ for constant 
respectively. The effect of the angular momentum barrier is more pronounced for small 77^ as in the case of the TH 
integrals with the regular Coulomb wave function as one may expect. The agreement between the approximations 
( [ZS] ) and ( [79| ) and the exact TH integrals shows the same trends as in the case of the TH integrals with the regular 
Coulomb wave functions. 
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FIG. 3: Scaled reduced Trojan-Horse integrals - 
of ^ = k/Q for angular mom enta Z = (a) and / 
approximations of Sec. [II C. 



1 (b). 



2 exp(7r77) for parameters t;^ = 0.01, . . . , 1.0 as a function 
The solid lines are the exact results and the dotted lines are the 



In the case of finite R the energy dependence of the TH integrals for f ^ can be extracted with the help of the 
approximations (|6^) and (^0|) which are vahd for 2ri ^ £_xq. The energy dependence of the differences 

jf^ (0, ry, fc, Q) - jf^ {R, r;, k, Q) = J^^^O, rj, - 4SA^o, V, (82) 

xo 

{2l + iy.Cii7^)^ ^da;v^/2i+i(2\/2^)^i(^) 



(2ry)'+i 
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0.1 



FIG. 4: Scaled reduced Trojan-Horse integrals Ji ,.^^{0,rj,^) exp{—-Kri) for parameters rj^ = 0.01, 1.0 as a function 
of ^ = k/Q for angular mom enta Z = (a) and I = 1 (b). The solid lines are the exact results and the dotted lines are the 
approximations of Sec. [II C. 



and 

j/^)(0, 77, k, Q) - J^^'^R, 77, k, Q) = j/fj,(0, 77, - 4%{^o,V, (83) 

XO 

2(277)'e 



(2Z + 1)!C,(77) 



is entirely given by the factor in front of the integrals since 77^ — rjk/Q is constant for constant Q. For large 77 and 
small ^ it is found that the energy dependence of the finite range correction to the TH integral with the regular and 
irregular Coulomb functions is given by ^^6*0(77) and £,/Co{r]), respectively. The finite range correction essentially leads 
to a change in the absolute value of the integral but the energy dependence of Ji'^^R, r/, k, Q) and Ji^\q, 77, k, Q) is 
the same for small k. 

When the Sommerfeld parameter rj increases, the coefficient Coij]) becomes very small and the TH integral with 
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Gi is much larger than the TH integral with Fi. The integrals are therefore dominated by Ji'^^ at low energies. 
In the case where a neutron is transfered in reaction the reduced TH integrals can be written as 



ih\ti{£.xo)ji{xo) - h\'^'{^xo)ji-iixo) 



in the variables xq — QR and ^ — k/Q with rj ~ 0. For ^ we find the approximation 

(±) (2/-1)!! . , , 

with a distinct /-dependence due to the centrifugal barrier. 



(84) 



(85) 



D. Simplified approximation of TH integrals 

The dominance of the contribution with the irregular Coulomb wave function G; in the TH integrals j/^'' at low 
energies also motivates a simple approximation in the limit fc ^ that has been used in applications of the THM 
to reactions with Ax ^ Cc. In this case the integral is completely determined by the contributions at radii close 
to R and the integrand decreases very fast with increasing r. If we assume that a general function /(r) decreases 
exponentially, the integral of /(r) from R to infinity can be calculated from 



I drJ{T)^f{R) j drexp(^~^)^ARfiR) 



(86) 



R R 

with AR = —f{R)/f'{R) and Rc(Ai?) > 0. In the TH integral j/^'' the appropriate value of AR is determined by 



AR. 



kG'i{r];kR) Qz[{Qr) 



Gii-q-kR) zi{Qr) 



(87) 



neglecting the small contribution of Fi for /c — > 0. From the approximation ([70|) of the irregular Coulomb wave 
function for kR <^ 2rj we find 



AR: 



2'qk K'^i+^{2^f^kR) Qz[{Qr) 
R K2i+i{2^2^) ziiQr) 



independent of k because ryfc is constant and Q is practically independent of k for small k. Applying this result to the 
TH integral leads to the approximation 



oo 

jl^\R,r],k,Q)^k J dr u[^\r];kr) zi{Qr) ^ kAR u[^\r];kR) zi{QR) 



(89) 



Since we are only interested in the fc-dependence, the TH integral can be replaced by the value of the integrand 

4'^\R,ri,k,Q) -.kRu\^\ri;kR) zi{QR) (90) 

at the cutoff radius R with the factor R instead of AR. The energy dependence for small k in this simple approximation 
agrees with the exact result for the TH integrals. 

IV. CROSS SECTIONS FOR SPINLESS PARTICLES 



In order to establish a closer connection of the three-body and two-body cross sections the Trojan-Horse T-matrix 
element (591) is recast in the form 



Tj/^{kcc,kBb;kAa) 



47r 



W{QBb) f"{QAa.kcc) 



kAxQAa 



(91) 
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with the Trojan-Horse scattering amphtude 



f^"iQAa,kcc) 



kAxQAa /^(-) 



An 



Cc, asym 



{kcc,rcc)\exp 



aA<Px 



(92) 



2ih 



^ V(2Z + l)Pl{kcc ■ QAa).l^ 
ikcc'^ V ^Ax 



S^xCcJl^'' (-R: VAx, kAx, Qao) - SaxCcJi ^ {R, VAx, kAx, Qao) 



which resembles an amphtude for the two-body reaction 

C + A + X . 



(93) 



The essential difference is the appearence of the Trojan- Horse integrals (^7|). The argument of the Legendre polynomial 
Pi in equation ( |9^ ) corresponds to the usual cosine of the scattering angle "dcm- It is convenient to decompose the TH 
amplitude 



J Jn + JC OAxCc 

into a nuclear part and a purely Coulomb part. The nuclear contribution 



rTH 

In 



2ikcc 



Y,{21 + l)Pl{kcc ■ QAa)^P^T'A^CcJi:^\R. VAx,kAx,QAa 



with 

'^AxCc = e 

depends on the nuclear S-matrix element 



VAx 



[SaIcc - SaxCc] 



qNI 

^a. 



xCc 



AxCc^ 



-icriivcc) 



(94) 



(95) 



(96) 



(97) 



which is obtained from the full S-matrix element by compensating the Coulomb phase shifts ai{r]cc) and (Ti{riAx) in 
the initial and final states, respectively. The Coulomb contribution 



rTH 

Jc 



^J2{21 + l)Pl{kcc ■ QAa)e'^'^'"'''^jf\R.riAx,kAx,QAa) 



(98) 



appears only if elastic two-body reactions are studied. It depends only in the TH integrals with the regular Coulomb 
wave function and is discussed in Section |v^ . 

If the Trojan-Horse integrals j/^'' are replaced by one the scattering amplitudes reduce to the standard form for 
the reaction (^3|). At low energies there are usually only few contributions to the nuclear scattering amplitude with 
small angular momenta I due to the increase of centrifugal barrier with I in the two-body reaction. 

The cross section ([l5| ) in the laboratory system can now be expressed in the form 



dEcd^cd^c 
with the two-body Trojan-Horse cross section 



W{QBb) 



TH 



2 167r^ vcc da 

^\xQ\a ""Ax dQ 



da 



TH 



dn 



'^Ax \jTH\^ 
VCc ' ' 



(99) 



(100) 



A corresponding expression holds for the cm. cross section ( |13| ) with the appropriate kinematical factor. This result 
has a similar structure as in the PWIA, i.e. a product of a kinematical factor, a momentum distribution and a two-body 
cross section. However, the momentum distribution \W\^ is not directly the ground state momentum distribution 
|$q|^ of the Trojan Horse a and the Trojan-Horse cross section ( |lO0| ) contains explicitly the off-shell effects. It is a 
sum 



da 



TH 



dn 



da^ 

dn 



dn 



SaxCc + 



da^ 

dn 



(101) 
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of a Coulomb contribution 



dVL 



— I fTH\ 

— \ JC 



a Coulomb-nuclear interference contribution 



daj 



TH 



dn 



E 



21 + 1 

kcc 



Pli'kcc ■ Qao) Re {lfE"y T'a.CcJi^\R, VAcc, kA^^QAa) 



and a nuclear contribution 



dn 



^ ^(2; + 1){21' + 1){10 I' 0|A 0)^PxCkcc ■ QAa) 



(102) 



(103) 



(104) 



T 



v* 



AxCc^ AxCc^l 



J, 



(+) 



{R 

7 TjAx J kAx J QAa J VAxi kj\xi Qao) 



Again, the expression for the TH cross section closely resembles the usual cm. cross section for the two-body reaction 
(p3|). The appearence of the TH integrals j/^"* accounts for the off-shell effects and the scalar product kcc ■ QAa 
appears as the argument of the Legendre polynomial instead of the cosine of the two-body cm. scattering angle "dcm- 
The TH Coulomb scattering amplitude (Bq) replaces the on-shell Coulomb scattering amplitude 



fc = 



1 



2ih 



Cc 



E 



(2/-|-l)Pi(cosi?c 



2k 



Cc sm 



2 



■ exp 



(exp[2icri(?7cc)] - 1) 
2icro(?7cc) - 2i77cclnsin 



(105) 



in the usual elastic two-body scattering amplitude. 

The expression for the TH cross section simplifies considerably in special cases. If Ax ^ Cc and only one partial 
wave I Ax contributes to the cross section the Trojan-Horse cross section is given by 



da 



TH 



dn 



dai 

'dn 



(Cc-> Ax) Pl{R,r]Ax,kAx,QAa) 



with the usual partial on-shell cross section 



dai 

dn 



{Cc Ax) 



(21 + 1)^ 



I 'T'^ I 
l-'-AxCc] 



A 



ilOlO\XO)^Pxikcc-QAa) 



for the two-body reaction (93) and the penetrability factor 

Pl{R,r]Ax,kAx,QAa) = Ji {R,VAx,kAx,QAa) ■ 

In the simple approximation of Sec. |lll D it is given by 

PiiR,Ti,k,Q)^eR^ [F^{ri;kR) + Cf{r,;kR)] zf{QR) 



(106) 



(107) 



(108) 



(109) 



From the fc-dependence of the TH integrals or the Coulomb wave functions a k^ exp{2TTri) behaviour of the penetrability 
factor is found for fc — > 0. For the transfer of a neutron the dependence of Pi on small k is given by fc^"^', see eq. 
(^), with an ^-dependence determined by the centrifugal barrier. 
With the theorem of detailed balance one obtains the direct relation 



d^<j 



dEcdncdnc 



lab 



W{QBb) 



167r^ vcc dai 
klcQ\a ^Ax dn 



{Ax Cc) Pl{R, riAx, kAx, QAa) 



(110) 



of the three-body cross section to the two-body cross section for partial wave This equation shows clearly the 
"parallelism" of the A-\-x reaction and the Trojan-Horse reaction: The cross section (IIC) is proportional to the cross 
section fo r the A + x reaction, modulated by the penetrability factor Pi. The factor Pi is directly related to the TH 
integrals ( |108D . It leads in general to an enhancement of the higher partial waves and it contains the exp(27r77) factor. 

A most convincing beautiful example of this parallelism is the comparison of neutron elastic scattering and the 
{d,p) reaction on -'^^N in the same energy range of the continuum in ^^N. There are the same peaks in both spectra, 
changed in magnitude according to the factor P; (ry = in the neutron transfer case). Also the s-wave resonance 
which appears as a destructive interference with the / = continuum shows up nicely in the (d, p) spectrum J15[ . 
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V. CROSS SECTIONS FOR PARTICLES WITH SPIN 



For the application of the THM to a particular reaction one has to consider that the nuclei participating in the 
reactions and (|^) usually carry a spin. In this case the definition of the relevant quantities becomes more intricate 
but the general procedure remains the same as in the spinless case. Spins of an individual particle i will be denoted 
in the following by Si with projection i>i. 

Here, the channel spin basis in the derivation of the expressions is used. In the two-particle reaction the spins sa 
and Sx in the initial state are coupled to the channel spin sax ■ Similarly, in the final state sc and Sc are coupled to the 
channel spin scc- Coupling scc with the spin of the spectator Sb gives the angular momentum s/ of the three-body 
final state. The initial state of reaction (j^) is characterized by the angular momentum Si which is obtained from 
coupling the spins sa and Sa- 

The cross section of reaction (||) in the laboratory system 

dEcdncdn. ^ {2SA + 1) E E E ^CcSf.f, kAa, S...)\' (111) 

SiS f I'il' f SCc 

is obtained by summing over final spin states and averaging over initial spin states with the three-body T-matrix 
element Tfi that carries now both momentum and spin indices. In the plane wave approximation it is given by 

Tf/^{kcc, qsb, sccSfVf, qAa, sm) (112) 
= (exp {iqBb ■ rsb) *i;5ym(Cc - 6, kcc, Sfi'f)\Vxb\ exp {iqAa ■ rAa) (j^AaisiVi)) 
with the asymptotic wave function 

*ym(Cc - 6, fccc, SCc, SfVf) (113) 
= E (''Cc ^Cc Sb l^blsf Vf)'i>is^^{Cc, kcc, SCcVCc)<Pb{sbVb) 



where 



■^\t}m{Cc,kcc.SCc,^Cc) (114) 

= #^EE E f C/jfl..,.(a;fc.^.)3^;"M°(a,r„)Z^-r---(fccc) 

as„ JM l^lc. " 

for Ta > R and 'I'iJym — for Tq, < R. The vector spherical harmonics 

y'fMia,f) =Y,il m s ''l-J MyYirn{f)Ms'^) (115) 
and the angular distribution functions 

ZjM{k) = Y.(^msi,\J M)Yi^{k) ■ (116) 

m 

are obtained by coupling the channel spins s with the corresponding orbital angular momenta I to the total angular 
momentum J. The radial wave functions 

and 

^tllc.sc. («; ko.ro.) = ^Sc..o. ^"''") (118) 

contain the general S-matrix elements Sf°f^f^ for a transition from a channel with quantum numbers IccScc in 
the partition Cc to the channel /qSq, in partition a. 
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The T-matrix in the plane wave approximation 

-* 4:71 -* -* ~* 

Tj"{kcc,qBb,SCcSfVf]qAa,Siyi) = W{QBb) ^'^"{QAa,kcc,SCcSfl'fSil'i) (119) 

again factorizes into a form with the momentum amphtude W and the scattering amplitude 

^^"{QAa,kcc,SCcSfl'fSiiyi) = (120) 
^IaTsaxIccSCc'^IaI^'^'' ^-^^1 ^Axi Qao) ~ ^AxCcSlAxlcc^SAxSCc'^lji^ i-^^VAxj kAxi Qao) 

which is a matrix in spin space. The angular distribution is determined by the function 

'^lixSAxlcoSCo (f^Cc, QAa, SfVfSiVi) (121) 
tEE E E E^'-™Cc(fccc)r,l™,JQ^,) 



47r 



2^ . - 



{Sx Vx Sb Vb\Sa l'a){SA VA Sq Va\si Vi){sCc VCc Sb V f) 

{I Ax rriAx SAx VAxU M){SA VA Sx yx\sAx VAx){lcc rtlCc SCc ^Cc\J M) 

= 47r ^ ^{sii^ij ni\sf i^f){lccmccj m\lAxmAx) 



with 



^u::2ic.scM^^^>''S^,SA) (122) 
= (-l)^--^^-^--^"^=L^^==y2i^V2i~+T 

/o„- 1 i\/o„ [ 1\ I Sa SAx I I SAx Sb Si II SAx J Ia 

(2, + l)(2s^. + l)| s.][ sf J sco]\lcc J sc.]- 
The Trojan-Horse scattering amplitude is a sum 

= ^n^+^c^^AxCc (123) 

of the pure Coulomb contribution 

•^C'' = E E E (2-/ + mLsAxlc.se. (kCc, QAa, SfUfSiUi) (124) 

J SAx lAxl'Cc 

e-^Ax irioc) jiF) ^^^^ f.^^^ QAa)SlAxlcJsAxSCc 

= /?^^.../..«,(-l)^»+^^+^''+^^^/2i^;+T^/2i:+T('/ ) 

1^ bf b}} ba J 

and the nuclear contribution 

= i;^^ ^ i'^J+mLsAxlc.saSkcc,QAa,SfiyfS„.,) (125) 

J SAx IaxIcc 
\ 'rr-^utsAxlCcSCc-^ilx (-^' ^^«) 

V ^Ax 

with 

rpJAxCc __ {riAx)+io-ir: Ancc) [qJAxCcN _ A < ^ A, , A 1 (^9a\ 

^lAxSAxlCaSCa - & Ax Cc I AxS Axl C.S c a AxC C^l AxI C a'^ S Ax SC .\ {^■^^) 

depending on the nuclear S-matrix element S/'^Z'^'^i^ » • 

^ ^ I^AxSAxI'CcSCc 
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The unpolarized cross section for the three-body reaction in the laboratory system 



dEcd^lcd^c 



{2SA + l)(2sa + 1) 

2 167r^ vcc do^ 



X! $Z X! l^/^C^Cc, ksb, SCcSfff, kAa, SiVi) 



(127) 



W{QBb) 



2 „,„ ^^TH 



is obtained by a summation over the spins in the final state and an averaging over the spins in the initial state. As 
in the case of spinless particles, the TH cross section 

dcr^^ 1 \ ^ \ ^ \ ^ I 

{2SA + l){2Sa + 1) ^,^.^^0^ 



-TH\ 



dQ 



(128) 



SiSf ViVf SCc 



dn 



daj 



TH 



dfl 



xCc 



da%" 

dn 



decomposes into a pure Coulomb contribution 



1 



dn 



{2SA + l)(2Sa + 1) 



nE EE I^^T = I/. 



TH\ 



iSf ViVf SCc 



an interference contribution 
daJ" _ 

dn 



{2SA + l)(2s„ + 1) 



nEEE2Re[^S--F; 



■THI 
N I 



(129) 



(130) 



SiS f ViVf SCc 



(2sA + l)(2sx + 1) 



lyEEE 



Re 



(*/< 



TlTsc:iccScAcl P^cc (kcc ■ Qau) 



Cc SCc 



2J+ 1 

kcc 



and a nuclear contribution 
da^ ^ 

dn 



1 \ ^ \ ^ \ ^ VAx I T-T/f 1 2 

{2SA + 1)(2S„ + 1) A. 2- 2^ VCc 



(131) 



SiSf ViVf SCc 



{2SA + 1){2S. + 1) 4k^ce ? § sitcc Sec Sec 



-iJAxCc 



rpj' AxCc* 



with factors 



'U^sj^JccSCc^V^^sa^I'^^sccJiVA^^ ^Ax, kAx, QAa)4JiR, VAx, kAx, Qao) 

Zi/ (As) = V(2J+1)(2J' + 1)(2; + 1)(2Z' + 1) (/ /' 0| A 0) I j, | 



(132) 



for the angular momentum coupling. The Coulomb contribution is the same as in the case of spinless particles. Notice 
the change from Sa to in the spin degeneracy factors. 

If only one orbital angular momentum I = Iax = I'ax contributes in the Ax partition to the inelastic two-body 
reaction A + x ^ C + c, the TH cross section 



da 



TH 



_ (2.Sc + l)(2.S, + l) rfcr; 

dn ^ {2s a + l)(2sx + 1) dn 
is directly related to the usual on-shell cross section 



{Cc^ Ax) Pl{R,'nAx,kAx,QAa) 



^(Cc^Ax) = ^ 
dn^ ' (2sc + l)(2se + l)4fc^, 



i.2 EE E E 



(133) 



(134) 



(_l)...-.c.^../ [\sAx)Ztiv^(\scc)Px{QAa ' ^Cc) 

rpJAxCc rpj' AxCc^ 
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with the penetrabihty factor (|108|). In (133) the only difference to the case with spinless particles is the appearance 
of the spin degeneracy factors. Applying the theorem of detailed balance the simple relation 



(Pa 



dEcdVlcd^lc 



W{Q 



Bb, 



167r^ vcc dai 



{Ax ^ Cc) Pl{R, rjAxi kAx, QAa) 



(135) 



is found which is the same as in the spinless case. The increase of the factor Pi at small kAx compensates the decrease 
of the two-body cross section ^ {Ax — + Cc) . 

VI. TROJAN-HORSE COULOMB SCATTERING AMPLITUDE 

In case of an inelastic two-body reaction only the nuclear contribution to the TH scattering amplitude is relevant. 

For the elastic scattering with Ax = Cc also the TH Coulomb scattering amplitude (98) contributes to the total 

( i^) 

TH scattering amplitude. It depends only on the TH integrals J; with the regular Coulomb wave functions. The 
summation over I in eq. ( p8[ ) poses no serious problem but it is convenient to reformulate the expression as in the 
case of the TH integrals in order to see the reduction of the Coulomb barrier more clearly. For that purpose the TH 
Coulomb scattering amplitude is written as a difference 



TH 



fTH _ f 
Jc — J CO 

of a contribution (suppressing the indices of the momenta) 



fTH 
JCR 



TH 

CO 



fc2 



{21 + \)Pi{k ■ Q)e'^'^'^'>M^\Fi,7^, k,Q) 



where the cutoff radius R is set to zero and a finite range contribution 



TH 
CR 



E 

I 



{2l + l)Pi{k-Q)e""'^'^^Li{R,r],k,Q) 



(136) 



(137) 



(138) 



with the integrals 



Li{R,r^,k,Q)^ dr Fi{-q-kr) zi{Qr) 



(139) 



These are easily calculated numerically and decrease rapidly with increasing / leading to a rapid convergence of the 
sum (138). Similar to relation ( |9^ ) for the full TH scattering amplitude, the contribution (137) to the TH Coulomb 
scattering amplitude can be expressed as a matrix element 



TH 



fco 



^(vI.W,(fc)|exp 



{iQ-f)) 



with the pure Coulomb scattering wave function 



(±) 

Coul 



(fc) — e 2T(1 ± ivi) e-Kp{ik ■ r)iFi{^iT], 1; ±i[kr ^k ■ r\) 



(140) 



(141) 



where li^i denotes a confiucnt hypergeometric function. This form allows an analytical calculation with the result 



fTH 
JCO 



27^PQ Co (7?) 

Q2 ~k^ {Q- kf 



exp < ictq ('?) + if] In 



Q^-k^ 

{Q~ky 



(142) 



as explained in Appendix |B[ Contrary to the ususal Coulomb scattering amplitude ( 105 ) there is no divergence for 
fc ^ since Q remains finite. In this limit, the energy dependence of the TH Coulomb scattering amplitude is given 
by fcCo(7?). The apperance of the Co{'r]) factor in the amplitude causes a strong reduction at small k. 



VII. APPLICATIONS OF THE TROJAN-HORSE METHOD 



Several reactions have been studied with the TH method recently. They are listed in Table | with ^H and ^Li 
(= a-|-d) as typical "Trojan Horses". These nuclei allow to study the transfer of protons, neutrons, deuterons and 
a-particles, which covers most of the cases of astrophysical interest for the two-body reaction. 
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TABLE I: Projectile energy -Epro and corresponding quasi-free energy E''J^ for pairs of two-body and Trojan-Horse reactions. 



two-body 
reaction 



Trojan-Horse 
reaction 



"^HfLi^H^ 



[MeV] 



[MeV] 



Ref. 



He 



^Li(p*He)*He 
i2cC'He,*He)i2c 
3He(2H,p)''He 
'^Li(p,^He)''He 



^Li( 

2 



■Li; 



He' 

4 



■He)' 

,4i 



He 



H(^Li,''He'*He)n 

®Li(=^He,p'*He)'*He 
2H(«Li,=*He-*He)n 



6.0 0.029 

19.0 - 21.0 0.161 - 0.412 

12.0 - 18.0 1.527 - 3.027 
5.8 0.845 
25.0 1.362 




A. Kinematical Conditions 

The Trojan Horses employed so far have a dominant s-wave configuration in their gound state. Their momentum 
distribution W{QBb) has a maximum around zero. Correspondingly, the equation 



iBb 



= (143) 



defines the so-called quasi-free condition. In this region of the three-body phase space the cross section for the 
quasi-free reaction will reach a maximum. From this condition the corresponding quasi-free cm. energy 

E'if^ = EAa(l-^^)~ea (144) 

in the initial channel of the two-body reaction (^ is derived from energy conservation (^l|) assuming the plane wave 
approximation. It is obvious that even with a large cm. energy EAa in the entrance channel of the three-body reaction 
(^ a small energy Eax can be reached if a suitable Trojan Horse a is chosen. This is confirmed in Table | where the 
projectile energy in the laboratory system and the quasi-free energy are shown for several reactions. 

The relation between i?^'^ and EAa is purely a kinematical consequence. It is not related to the Fermi motion of 
particle x inside the Trojan Horse a which would involve a dependence on the width of the momentum amplitude. 
However, in an actual experiment a cutoff in the momentum transfer Qsb is chosen to select the region where the 
quasi-free process dominates the cross section. This procedure limits the range in energies Eax that are within reach 
in the experiment for a chosen projectile energy -Epro- 

In case of the quasi-free condition, all nuclei in the final state of reaction (^ are emitted in the same plane. The 
momentum of the spectator h is in beam direction which makes it difficult to detect b in the experiment. In the 
laboratory system the nuclei C and c are emitted under angles dc and on opposite sides of the beam axis. If the 
scattering angle ^cm in the cm. system of the two-body reaction is given, the angles t9c, the so-called quasi-free 
angles, are completely specified for a fixed beam energy from kinematical considerations. 

As a consequence, the quasi-free condition determines the setup of the experiment. If a particular two-body reaction 



(g) is to be studied close to a cm. energy Eax and if the Trojan Horse a is selected, then from equation (144) the 
appropriate beam energy can be extracted. Since the momentum amplitude W has a finite width it is possible to 
study the two-body reaction in a certain energy window around Eax- The cm. scattering angle "dcm determines the 
arrangement of the detectors close the pair of quasi-free angles. 
When EBb is zero in (^) the maximum energy 

E^^'' = EAa - Sa (145) 

in the two-body system is reached for a fixed cm. energy Eau in the entrance channel of the three-body reaction. In 
this case ksb = and QAa = kAa- The momenta of all nuclei C, c, and b in the final state are parallel to the beam 
momentum in the laboratory system. Since £a > it follows that the relation 



kAx<J — QAa<QAa (146) 
fJ-Aa 



holds in all kinematically allowed regions of the phase space. 
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B. Threshold behaviour of cross sections 



The energy dependence of the two-body cross section 



%iA.^Cc)= (147) 



^— E E(2^ + + ^' 0| A 0)2p, (cos ^)^Lcc^a:cc 



for the inelastic A + x ^ C + c reaction above the reaction threshold is governed by the 1 / factor and energy 
dependence 

SaxCc « exp{-Trr]Ax) (148) 
of the relevant S-matrix element. The Coulomb barrier leads to a strong suppression of the the cross section 

^(Ax ^ Cc) oc k^l exp{-2TTr]Ax) (149) 

for kAx — >■ due to the decreasing exponential factor. This behaviour motivates the introduction of the astrophysical 
S-factor (|^). In the TH cross section that appears in eq. (^), the factor k^l is replaced with k^^ and the TH in- 

tegrals j/^^ appear. Their energy dependence for small kAx is determined by kAx/Co{r]Ax) ~ kAx GxpiT^VAx) / \/^t^VAx 
from the contribution of the irregular Coulomb wave function. This leads to a kAx dependence of the three-body 
cross section (^9|) according to 

dEctncdn^ « kfxVAl ^-P(-2-^A.)fcL^^|^^^ = i2nrjAxVAxr' = const. (150) 



in the lowest order of kAx', cf. also eq. (110) with the kAx dependence of the penetrability factor Pi oc fc^^ e'x.p(2Trr]Ax)- 
As a result the cross section does not vanish at the treshold but takes on a finite value. Of course, the same 
considerations apply to the case when the spins of the particles are considered. Also in the case of the transfer of 
neutron, like in a (d,p) stripping reaction, it is well known that the cross section is finite at the threshold En = 
[ pj| , |l8|. The reason is the same as in the case of charged particles. 

C. Extraction of astrophysical S-factors 

In an actual TH experiment the measured cross section of the three-body reaction depends on the geometry of the 
s etup , the chosen Trojan Horse, and the energy of the projectile. The differential three-body cross section ( |99| ) or 
( |llO| ) can be projected onto a simple cross section 



depending on the A + x cm. energy E. The efficiency function F takes cut-offs in particle energies, momenta (e.g. 
Qsh), the detector geometry etc. into account. The experimental spectrum ( |l5l| ) can be compared to a corresponding 
theoretical cross section from a simulation of the experiment assuming a certain energy dependence of the relevant 
on-shell S-matrix elements of the two-body reaction, e.g. from a R-matrix parametrization. By a variation of the 
parameters, the best fit to the experiment is found and the on-shell two-body cross section can be calculated. If there 
i s on ly a contribution of one partial wave, the procedure becomes simpler. The ratio of the measured cross section 



(151) to the corresponding simulated cross section directly gives the energy dependence of the S-factor relative to 
the energy dependence assumed in the theoretical S-factor of the two-body reaction. Due to the uncertainties in the 
description of the reaction mechanism it is expected to be difficult to extract absolute values of the cross section. 
So it is better to normalize the cross section to results from direct experiments at higher energies. However, we 
expect that the energy dependence of the cross section (S-factor) can be extracted much more reliably. For recent 
applications of the TH method with detailed information on the experimental realisation we refer to the references 
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D. Elastic scattering with the Trojan-Horse method 

The main aim for applying the TH method is the extraction of the energy dependence of cross sections or astro- 
physical S-factors for inelastic two-body reactions. But also the indirect investigation of elasting two-body scattering 
A + X A + X can be rewarding. In the direct t wo-b ody scattering process the cross section is dominated by the 
contribution of the Coulomb scattering amplitude ( |105| ) at low energies that diverges with fc^^. By way of contrast, 
the TH Coulomb scattering amplitude (^) vanishes for k^^ due to the appearance of th e TH integrals with 



the regular Coulomb wave function Fi. The nuclear contribution to the TH cross section (101) becomes dominant 



because of the TH integrals with the irregular Coulomb wave functions in eq. (104). This allows the study of nuclear 
effects in the scattering at small energies, e.g. the influence of sub-threshold or low energy resonances. First attempt 
in this direction were made in recent experiments p7| , 

VIII. SUMMARY AND OUTLOOK 

In this paper the basic theory of the Trojan-Horse method was developed starting from a distorted wave Born 
approximation of the T-matrix element. The essential surface approximation allows to find the relation between the 
cross section of the three-body reaction and the S-matrix elements of the astrophysically relevant two-body reaction. 
In the modified plane wave approximation the relation between the three-body and two-body cross sections becomes 
very transparent. The three-body cross section is a product of a kinematical factor, a momentum distribution and 
an off-shell two-body cross section. Off-shell effects are expressed in terms of so-called TH integrals that were studied 
in detail. Their energy dependence leads to a finite cross section of the three-body reaction at the threshold of the 
two-body reaction without the suppression by the Coulomb barrier. This allows to extract the energy dependence of 
astrophysical cross sections from the three-body breakup reaction to very low energies without the problems of electron 
screening and extremely low cross section. A comparison of results for S-factors from direct and indirect experiments 
can improve the information on the electron screening effect. However, dedicated Trojan-Horse experiments are 
necessary in order to achieve a precision comparable to direct measurements. 

The validity of the Trojan-Horse method can be tested by comparing the cross sections extracted from the indirect 
experiment with results from direct measurements of well studied reactions. In principle it is possible to assess system- 
atic uncertainties of the Trojan-Horse method by studying various combinations of projectile energies, spectators in 
the Trojan Horse and scattering angles. Furthermore, different theoretical approximations can be compared, e.g. full 
DWBA calculations with and without the surface approximation and simpler modified plane wave approximations. 

One may also envisage applications of the Trojan-Horse method to exotic nuclear beams. An unstable projectile hits 
a Trojan Horse target allowing to study specific reactions on exotic nuclei. A study of low-energy elastic scattering 
with the Trojan-Horse method opens another application which can lead to improved information relevant to the 
theoretical description of nuclear reactions at low energies. 
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APPENDIX A: ANALYTICAL CALCULATION OF RADIAL INTEGRALS 

The radial integrals 

oo 

Mu\H,7^,k,Q) = drHiir,:kr)r-'ziiQr) (Al) 





appearing in equation (|5^) can be calculated explicitly by using the integral representation 

i(fcr)'+i 

{21 + ly.Ciifi) 



^ ^ i-\-oo 

Fi{T^-kr)+iGi{T]-kr) = , , / e-'='-^(s - 2)'-"'(s + (A2) 
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of the combined Coulomb function Hi = Fi+ iGi . This form has been obtained from the representation in |Q by the 
substitution t — kr{s — i). The integration path in the variable s is a parallel to the real axis in the complex plane. 
Since the recursion relation (^) connects integrals of three successive Z-values, only the basic integrals with / = 
and I — 1 are needed for the explicit calculation of MjJ^ for all I. With zo{x) = sinx and zi{x) = sinx/x — cosx the 
integration over the radial coordinate r is easily performed and one finds 



i+oo 



Moo\H,7j,k,Q) = 



Coin) 



(s + i)"'(s-i)-*'' 

fc2s2 + Q2 



and 



M[-,\H,r,,k,Q) 



ikQ 



2+CO 



3Ci{r,) 

The substitution s = i{t -\- \)/{t — 1) leads to the expressions 

M-\H,T^,k,Q) = 

M-,\H,T^,k,Q) = 

with the integral 



ds 



1 — ir} 



-2 



d 



Co(r/)(Q2_fc2 
AkQ 

3Ci(??)(g2 - fc2)2^ 
h{-n,z) ^ I dt 



t2 ~2tz + l 



(A3) 

(A4) 

(A5) 
(A6) 

(A7) 



The contour C is the straight line from 1 to 1 + zoo parallel to the imaginary axis in the complex plane. The integral 
depends on the variable 



(A8) 



which is always larger than 1 for the conditions of the TH method. The integrand has two poles at zi = z+\/ z"^ — 1 > 1 
and Z2 — z — \J z^ — \ < 1 on the real axis. In the next step the path of integration is deformed to lie on the real axis 
and the integrand is broken into partial fractions. This yields 



dt 



1 



— ; — m Res2, — 

<2 ~2tz + l ' i2 _ 2tz + 1 



(A9) 



Zl - Z2 

with the contribution of the residue at the pole zi. The remaining integrals 



L = / dtt"^ 



1 



t - Zi 



for j = 1, 2 are evaluated with the help of the geometric series and the relation 



Trr) 



tanh(7r7]) 



oo 



n—l ' 



,2 ■ 



Collecting all contributions we find 



TTT] 



(AlO) 



(All) 



(A12) 



for the fundamental integral ( [AT] ). Separating real and i mag inary parts, the explicit expressions of the integrals 
M^q{Hq^ 7], k, Q) and M^f^{GQ^ri, k, Q) are found with eq. (A5). In the case / = 1 the derivative of Io{z) with respect 
to z has to be calculated first. 
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APPENDIX B: TROJAN-HORSE COULOMB SCATTERING AMPLITUDE 



Recalling the procedure to derive the radial integrals (Al), the contribution to the Coulomb scattering ampli- 
tude can be written as 



rTH _ 

27r(fc2 - Q 



2y(*cli(^)k"Vxp(4-rj) 



(Bl) 



by considering the Schrodinger equations for the Coulomb scattering wave and the plane wave. The matrix 

element can be evaluated by employin g the integral representation of the confluent hypergeometric function in 
the Coulomb scattering wave function (|l4ll). After the integration over the spatial coordinates, an integral remains 
that represents a hypergeometric function [[34| . This technique is similar to the evaluation of Bremsstrahlung matrix 
elements |42, |3|. One obtains 

(*cli(fc)k"'|cxp (*g-r)) = 

with the argument 

X = -2 



47re 
k- (Q-k) 



r(l + 277)2^1(1, -ir;;l;x) 



(B2) 



(Q - fc)2 

in the hypergeometric fmiction 2F1 which reduces to the simple form 

2Fi(l,-z77;l;x) = (l-x)'" 



(B3) 



(B4) 



for the given parameters. Combining the above results the scattering amplitude assumes the form (142) when the 
relation 



r(l + if]) = Co{ri) exp (^^rj + zero) 



(B5) 



with the Coulomb phase ctq is used. 

One application of this formula is the calculation of the modified momentum amplitude ( ^) with a Coulomb 
scattering wave function. In this case one finds 



W{q,k) = -M/(<z)^^^|^exp|7ao(^)+z,;ln 
where the indices Bb of rj and k have been suppressed for clarity. 



(g-fc)2 



(B6) 
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